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Abstract 

In this paper we investigate a class of solutions of Einstein equations for 
the plane- symmetric perfect fluid case with shear and vanishing accelera- 
tion. If these solutions have shear, they must necessarily be non-static. We 
examine the integrable cases of the field equations systematically. Among 
the cases with shear we find three classes of solutions. 

KEY WORDS : Exact solutions, plane symmetric models, inhomogeneous uni- 
verse 

PACS No. : 04.20.-q 

1 Introduction 

The standard Friedmann-Robertson- Walker (FRW) cosmological model pre- 
scribes a homogeneous and an isotropic distribution for its matter in the de- 
scription of the present state of the universe. At the present state of evolution, 
the universe is spherically symmetric and the matter distribution in it is on the 
whole isotropic and homogeneous. But in its early stages of evolution, it could 
have not had such a smoothed out picture. Close to the big bang singularity, 
neither the assumption of spherically symmetric nor of isotropy can be strictly 
valid. From this point of view many authors consider plane symmetry, which is 
less restrictive than spherical symmetry and provides an avenue to study early 
days inhomogencities. In recent years, cosmological models exhibiting plane 
symmetry have been studied by various authors (Rcndall, 1995; Taruya and 
Nambu, 1996; Zhuravlev et at, 1997; da Silva and Wang, 1998; Ori, 1998; An- 
guige, 2000; Chervon and Shabalkin, 2000; Nouri-Zonoz and Tavanfar, 2001; 
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Yazadjiev, 2003; Pradhan et at, 2003, 2005; Saha and Shikin, 2004, 2005) in 
different context, fnhomogeneous cosmological models play an important role 
in understanding some essential features of the universe such as the formation 
of galaxies during the early stages of evolution and process of homogenization. 
The early attempts at the construction of such models have done by Tolman, 
(1934); and Bondi (1947) who considered spherically symmetric models. Inho- 
mogeneous plane-symmetric model was first considered by Taub (1951, 1956) 
and later by Tomimura, (1978); Szekeres, (1975); Collins and Szafron, (1979a, 
1979b); Szafron and Collins, (1979). Bali and Tyagi, (1990) obtained a plane- 
symmetric inhomogeneous cosmological models of perfect fluid distribution with 
electro-magnetic field. 

In this paper we study plane-symmetric inhomogeneous models in presence 
of perfect fluid with shear and vanishing acceleration. Among the cases with 
shear we find three classes of solutions. 



2 Field Equations 

We consider the plane-symmetric line element in the general form 

ds 2 = -A 2 (r,t)dt 2 + B 2 (r,t)dr 2 + C 2 {r,t)(dx 2 + dy 2 ). (1) 
The four velocity of the fluid has the form 



0,0,0,1 ) 



(2) 



The energy-momentum tensor in the presence of a perfect fluid has the form 

Tab = (p + p)u a ii b + pg ab , (3) 

where p and p are the energy density and the pressure of the fluid, respectively. 
In this coordinate system the Einstein's field equation 
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The dot denotes partial derivative with respect to time, the prime indicates par- 
tial derivative with respect to the coordinates r. k is the gravitational constant. 
The consequences of the energy momentum conservation 
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are the relations 
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The plane symmetric solutions can be classified according to their four kinemat- 
ical properties, i.e. rotation, acceleration, expansion and shear. In the comoving 
frame of reference these quantities read 



Uab = U[a-b] + U[o"b] = 0, 
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The square bracket in Eq. pip denote antisymmetrization and the round 
bracket in Eq. fTJJ denote symmetrization. 



3 Plane symmetric Non-static solutions with shear 
and without acceleration 

To simplify the field equations we assume vanishing acceleration of the four 
velocity. According to Eqs. (|10[) and (fl2|) the pressure of the fluid and the 
metric function In (A) depend on time only, i. e. 

P = P{t), A = A(t) (16) 
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Eq. 



reads now 



& 
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We have to investigate two essential different cases namely 

C = 0, 

and the general case 

C = C{r,i). 

First of all we investigate the case C = C(t). 



(17) 

(18) 
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3.1 Solution for the case C = C(t), B = B(t) 

Without loss of generality, we can take c = t. The field equations take the form 



Kp = T 2 



1 2B 
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Eq. (|20[) defines the energy density. Eqs. (|21|) and (|2"2"|) define the pressure. The 
only equation which actually has to be solved is the condition of consistency of 
Eqs. ([SI]) and ||22J) 

f(t)v = g(t)v + 1. (23) 
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Here we have used the following abbreviations: 
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, , , Bt Bt 2 
ff(*) = l-B— 
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Eq. (f25|) contains the two unknown quantities A = A(t) and B = B(t). We get 
all solutions of this class by choosing the function B = B(t) arbitrarily. Then 
we determine the function v = by solving the linear first order Eq. ([23]) . 
The solution reads 
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The energy density, the pressure, and the shear can then be written as 
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Out of the infinite number of exact solutions of Kantowski-Sachs class we present 
only the solution following Mc Vittie and Wiltshire (1975). Let us choose 

B = n\nt, (31) 

where n is constant, n 2 ^ 1. The functions f(t), g(t) and v(t) read 



1 



/(f) = -t(n-l), 
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where Co is an arbitrary constant. In this case the geometry of the universe (fTJ) 
takes the form as 



ds 2 



(n 2 



1) 



dt 2 + {n\ntydr 2 +t 2 (dx 2 + dy 2 ). (35) 



_l + C*o(n 2 -l)i- 2 ( n+1 )_ 
The energy density, the pressure and the shear for the model (I35|) are given by 
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From Eqs. <PS|) and (|3T|) . we observe that p > 0, p > and p is a decreasing 
function of time. This model satisfies the weak and strong energy conditions 
and also has a physically acceptable fall-off behaviour in both r and t. 
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3.2 Solution for the case C = C(t), B = B(r,t) 

In this case we write down the Einstein's field equations as follows: 
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From Eqs. (|40|) and (|41[). by using the condition of consistency, we have 
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Equation (|42p is a non-linear partial differential equation. To solve this equation, 
following the technique of Hcrlt (1996), we choose 
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where Co is an integrating constant. Using Eq 

B + PB - NB = 0, 

where 



we obtain 
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Eq. (|44|) has the special solution B = t n . Following D'Alcmbcrt's method. We 
insert the expression 



B = D(r,t)t n 

in to Eq. (|44j) . where D(r, t) is a function of r and t, wc obtain 
D+{P + 2nr 1 )D + (n - 1)P(1 + nt^D = 
Integrating Eq. (|48]) . we obtain 
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+ C 2 e--t- n WhittakerW 
From Eqs. (@7J) and (gU), we have 
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The energy density, the pressure, and the shear can be computed easily now by 
means of formulae (O, © and (fTBj) . The geometry of the universe ((T])in this 
case reduces to 



ds 2 = -A 2 dt 2 + B 2 dr 2 + t 2 {dx 2 + dy 2 ), 
where A and B are given by (|4*3")l and ([50]) 



(51) 



3.3 Solution for the general case C = C(r,t), B = B(r,t) 

The pressure p and the function A depend on the time t only. By means of a 
scaling of the coordinate t we arrive at 



A = l, C 



BC 



In this case we write down the field equations as follows: 
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The second equation for the pressure follows after differentiation of Eq. (fB"4"|) . 
We can integrate Eq. (|55|) to obtain 



In C = 1115 + -In [l-eh 2 {r)\. 



(56) 



Here | In [1 — eh 2 (r)] is the integration "constant", where %{r) is an arbitrary 
function of r. We get 



(57) 
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The remaining equation to be solved is ([M]) . We use the result ([57)) and get 

p(t) = ^[l-eh 2 (r)-{2CC + C 2 }] : (58) 
with the abbreviation 

C(r,t) = zi(r,t), (59) 



we obtain from Eq. ([58]) the nonlinear differential equation 

Z +^k P Z =^(l-eh 2 (r)) Z~i. (60) 
The mass density follows from eq. ([53]) by the use of ([57]) and ([58]) 

k P = ^ [(1 - eh 2 (r)) C" + h{r)h'{r)C] - k V (2C"C + C' 2 )+ 
• 2 /2C" C' 2 1\ 2C'C 2C' 2 \ 

c {— + ^ + ^j + w +c ( 4C + —) (61) 

To obtain solutions we have to proceed as follows. First of all we choose some 
p(t) ad hoc. Then we have to solve Eq. (|60[) . The last step is to the calculation of 
mass density and shear. Eq. (|60[) can be simplified by means of the abbreviations 

3 / 1.^2, 



Z(r,t) = EZ{r,t), E* = - I -- + ft» ) . 
Then we get 

Z+ -npZ = -eZ~z (63) 

with 



C(r,t) = ^(~+ft a (r))^- (64) 

The regularity conditions at the base of the plane are fulfilled automatically if 
the arbitrary function h~(r) is chosen as 

k(r) =0 at r = 0. (65) 

In this case Eq. ([6"!) reduces to 

*P = ^ - -P(2C"C + C' 2 ) + C 2 + + i_ 

2C'C . ^/.^„ . 2C 12 



C'C \ C 



C AC" + — (66) 



At first glance the solution of ([63]) seems rather trivial. We can choose an 
arbitrary function Z = Z(t) and then determine the pressure. But this does not 
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work. First of all notice that in the case Z = Z(t)a(r) the function C(r, t) has 
the form 



C{r,t) = J\ (-\ + h\r) 



As a result, the shear vanishes because 



1 2 1 3 1 C B\ Id, 

a 1 = a > = V , = 3 c-5r3i ,ln 



a_ 
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(67) 



(68) 



Therefore Z has to depend on the co-ordinate r in the form 
Z = Z{Ji(r),J 2 (r),i} fore^O. 



(69) 



Here the arbitrary functions J\ (r) , and J2 (r) are integration "constant" which 
enter the general solution of Eq. ([63]) . We see that the r-dependence in the 
solution comes from these integration "constant" . 
If we set e = — 1 and %(r) = 0, then by (|56[) 



C = B. 

Hence, in this case, the metric ([I} reduces to 



ds 2 = -dt 2 + C' 2 (r,t)dr 2 + C 2 (r,t)(dx 2 + dy 2 ). 



(70) 



(71) 



For model (|7T|) , the corresponding expressions for pressure and density can be 
computed from Eqs. (|63|) and ((66)) respectively. 



Let us recall here two important families of solutions corresponding to par- 
ticular choices of p(t) (we do not claim to be exhaustive). 



3.3.1 The case p = constant 7^ 

In this case the differential equation (|63[) reduces to the form 

Z + NZ = -eZ~§ (72) 
where N = \np. Integrating ([72]) . we get expression for Z 

N ~ 9-5 

jZ 3 + — eZ-z +(l-k 1 )Z-k 2 = 0, (73) 

where k\ and k% are integrating constant. Putting the value of Z from (|73p in 
(|6"T|) . we get the value of C(r, t) which is not reported here due to complexity 
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3.3.2 Dust case p — 

In this case the differential equation (|63p reduces to the form 

Z = -eZ~^. (74) 
Integrating (|74|) . we get expression for Z 

(1 - mi )Z + ^eZi -m 2 = (75) 



where mi and m2 are integrating constants. Putting the value of Z from (|75[) 
in ([57]) . we get the value of C(r, i). 

3.3.3 The Explicit Form of the General Solution 

As stated in the preceding section, a solution of field equations (|53|) - ([55"1) can 
be obtained for any given form of p(t). We give here the explicit form of the 
general solution with no restriction on p(t). 

Let us perform in (|63[) the following substitution 

e = (76) 

so that we obtain the solution 

Z + -^npZ = 0, (77) 



which is linear in Z . The general equation of (|77|) is of the form 

Z = D ia (t) + D 2 b(t) for e = 0, (78) 

where Di and D 2 are arbitrary functions of r and a and b are two different 
(nontrivial) particular solutions of (|77[) . In fact, we need only know one non- 
trivial particular solution of {77]), namely a(t). This is because both a and b 
verify (|77j) . that is 

3 

a+ -np(t)a = 0, (79) 

and the analogous expression for b(t). Both functions must then be related to 
one and another by 

-a - I (80) 
a b 

which can be solved, up to a quadrature, for b once a(t) is known. A particular 
choice of b(t), different from a(t), is given by 

b(t) = a(t) J a- 2 dt' (81) 

Note that ([79]) actually may be interpreted as a mere definition of p(t) (because 
a depends only on time). We can consider either p(t) being arbitrary and a(t) 
derived through (|79[) . or vice verse. 
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4 Conclusion 



We have studied a new class of plane symmetric inhomogcneous cosmological 
models with shear and vanishing acceleration in presence of a perfect fluid. Our 
analysis is an attempt to obtain more exact solutions so that our understand- 
ing of these objects may be improved. It is hoped that some of the solutions 
presented here will provide the basis for a detailed physical analysis of plane 
symmetric inhomogeneous models with shear and vanishing acceleration in the 
gravitational context. The solutions we had obtained for the case p ^ and 
p = (dust case) and explore the consistency of the formalism and outlined 
obtaining generalization. Some of the models are being persued at present and 
would be reported in future work. 
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